Transverse force on a quantized vortex in a superconductor 
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The total transverse force acting on a quantized vortex in a type-II superconductor determines 
the Hall response in the mixed state, yet a consensus as to its correct form is still lacking. In this 
paper we present an essentially exact expression for this force, valid in the superclean limit, which 
was obtained by generalizing the recent work by Thouless, Ao, and Niu [D. J. Thouless, P. Ao, and 
Q. Niu, Phys. Rev. Lett. 76, 3758 (1996)] on the Magnus force in a neutral superfluid. We find 
the transverse force per unit length to be f = p K x V, where p = p n + p B is the sum of the mass 
densities of the normal and superconducting components, K is a vector parallel to the line vortex 
with a magnitude equal to the quantized circulation, and V is the vortex velocity. 
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The discovery of high-temperature superconductivity 
has stimulated a renewed interest in the statistical me- 
chanics and dynamics of vortices in type-II superconduc- 
tors [jl[. Despite the tremendous progress made, however, 
the answer to one of the simplest questions, namely, the 
form of the equation of motion for a single isolated vortex, 
remains controversial. If we let R.(t) denote the position 
in the xy plane of the center of a straight vortex line as 
a function of time, then the classical equation of motion 
can be written as 



M 



= -v 



rfR 
~dt 
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~dt 



xe, + f p (R) + f d (R). (1) 



Here we have taken the circulation vector K of the vortex 
(a vector parallel to the vortex with a magnitude equal to 
the circulation Q|) to be along the z direction. The first 
two terms on the right-hand-side of (]!]) are to include all 
forces linear in the vortex velocity. Here f p represents the 
pinning force due to disorder (averaged over an area set 
by the size of the vortex core), and fd includes the vari- 
ous "driving" forces possible, such as the Lorentz force, 
which may depend on the density of the normal and su- 
perconducting components n s and n n , and therefore on 
the vortex position R, but do not depend on the vor- 
tex velocity. The equation of motion (Q) also describes 
a vortex in a Bose or Fermi superfluid, where f p might 
describe the force from an externally imposed wire as 
in a Vinen-like experiment ||] , and fd would include the 
superfluid-velocity-dependent part of the Magnus force 
^ and possibly other vortex-velocity-independent con- 
tributions §. 

The difficulty concerns the determination of the co- 
efficients M, T], and 7, which describe the vortex ef- 
fective mass per unit length, viscous damping force per 
unit length, and nondissipative transverse force per unit 
length, respectively. Of particular interest and especially 
controversial is the coefficient 7, which in a neutral super- 
fluid has been recently shown by Thouless, Ao, and Niu 
H (hereafter referred to as TAN) to be 7 = p s K, where 



p s is the mass density of the superfluid component far 
from the vortex core, and K is the quantized circulation. 
For a vortex in a neutral Bose superfluid with q quanta 
of circulation, this coefficient is simply 7 = qhn s , with 
n s the superfluid number density away from the core, 
whereas in a (paired) Fermi superfluid 7 = ^qhn s . The 
TAN result shows that the transverse force is universal, 
independent of the detailed microscopic structure of the 
vortex and its interaction with the normal component of 
the fluid. The universality of the nondissipative trans- 
verse force is not found in some approximate microscopic 
calculations, however, which obtain coefficients that re- 
duce to the TAN result only in certain limits 

In this paper we calculate the nondissipative trans- 
verse force on a quantized vortex in a charged superfluid 
or superconductor by following the method of TAN ||. 
Initially, we shall consider a system with translational 
invariance, and to address a more realistic model of a su- 
perconductor we include afterwards a periodic potential 
from the lattice, assuming the lattice constant b is small 
compared with the coherence length £. In both cases we 
find the transverse force per length to be given by 



f = pKxV, 



(2) 



where p — p n + p s is the total mass density of the fluid 
away from the core, K is the quantized circulation vec- 
tor defined above, and V is the constant vortex velocity. 
Most relevant is the fact that, as in TAN, our results for 
this superclean limit of a superconductor also indicate 
a universal vortex-velocity-dependent part of the trans- 
verse force. In the notation of Eqn. (|]), we find 7 = pK. 
For a vortex with q trapped flux quanta $0 = hc/2e we 
can write this as 7 = ^qhn, where n is the total density 
of the fluid far from the core. The reason that the trans- 
verse force in the charged case differs from that for the 
neutral superfluid, which has a coefficient p s instead of 
p, is related to the Meissner effect, as will be explained 
below. 
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The Hamiltonian for the charged superfluid or super- 
conductor is taken to be 



where 



(3) 



(4) 



is the Hamiltonian for N bosons or fermions with neg- 
ative charge — e, interacting with each other and with 
a uniform positive background charge, and (in Gaussian 
units) 



2to 2 ( 



£K^'(r„-r„0^, 



(5) 



with T ij (r) = (S'^lr^ 1 + rV|r|~ 3 )/2c, is the current- 
current interaction, which, as first noted by Darwin [ pl)[ , 
correctly accounts for the electrodynamics in the trans- 
verse gauge through order v 2 /c 2 . The U in (Q) is a 
Coulomb interaction term, and in the case of a super- 
conductor also contains an additional short-ranged at- 
tractive interaction to produce superconductivity. Also, 
Hb = —J2nf e 2 h\r n — r| _1 accounts for the inter- 
action with the uniform positive background charge en. 
The full Hamiltonian also contains a cylindrically sym- 
metric pinning potential V centered at position R p| . 

Before proceeding it is important to note that the 
Hamiltonian (||), with the current-current interaction 
term, provides an accurate model of a charged super- 
fluid or superconductor. In particular, we show in the 
appendix that it exhibits a Meissner effect with the con- 
ventional London screening length. 

We shall follow TAN and calculate the force on the 
vortex by expanding the time-dependent wave func- 
tion, given by ihd t \V(t)} = H(R(t))\V(t)), in a basis 
of instantaneous eigenstates satisfying H(R) \ip a (R)) = 
£ a (R)|V> Q (R)). We shall assume a three-dimensional 
system with the z direction along K. The system is 
assumed to be infinite in the x and y directions, and 
L is the thickness in the z direction, which is also the 
length of the vortex. In the absence of the pinning 
potential the system is therefore translationally invari- 
ant in the x and y directions. The translational in- 
variance allows the instantaneous eigenfunctions to be 
taken as ip^(ri, • • • , r/v) = ip a (^i — "■, fjv — R), where 
tp a (ri, ■ ■ ■ , r^r) arc the eigenfunctions with the pinning 
potential centered at the origin. Initially, at a time to, a 
vortex is assumed to be bound to the pinning potential 
at R(to) = Ro m a state characterized by the density 
matrix J2 a fa\ipa(^o))(ipa{^o)\, where f a is the occu- 
pation probability of state a. The force on the vortex 
is F = -£ a /„ {* a (t)\V R H\* a (t)), where |* (t)> is 
the solution of the time-dependent Schrodinger equation 
starting out in \ip a (R.(to))) . Following Ref. || we obtain 
a transverse force equal to 



^(e z xV) jU.(V-V')p(r,r')L 



(6) 



where p(r,r') is the one-particle density matrix and L is 
the thickness of the system in the z direction. In (|^) we 
choose the radius a of the integration contour to be much 
larger than the London penetration depth A. 

The integrand in (^|) is proportional to the canoni- 
cal momentum density. The current-current interaction 
term (|5|) has played no role up to this point, and, indeed, 
the expression (|6|) in the two-dimensional limit is identi- 
cal to that obtained by TAN for the neutral superfluid. 
However, the resulting force is different, as can be seen by 
writing (^) in terms of the canonical momentum density, 
which leads to a transverse force per unit length equal to 



f = e 2 x V 



/dl-j(r). 

J a 



(7) 



The TAN result for neutral superfluids is obtained by 
writing the momentum density as j = p n v n + p s v s and 
assuming that the viscous normal component does not 
circulate |12|. In the case of a charged superfluid or su- 
perconductor, however, the gauge-invariant momentum 
density j + -nA vanishes in the region containing the 
integration contour because of the Meissner effect (see 
Appendix), and we therefore obtain (^). 

Our conclusion suggests that both the normal and su- 
perconducting components of the fluid contribute to the 
total transverse force. However, we would like to empha- 
size that the model considered here does not include dis- 
order, and, as such, leads to a normal component having 
an infinite conductivity. This neglect of disorder may be 
responsible for the p dependence of the transverse force 
instead of the usual p s . 

Finally, we would like to remark that the result (||) 
holds even in the presence of a lattice, as long as the 
coherence length or vortex core size £ is large compared 
with the lattice constant b. To demonstrate this we add a 
periodic potential term v(r n ) to Hq and expand the 
instantaneous eigenstates in a basis of localized Wannier 
functions a(r — 1) for the relevant band, 

^(r 1 ,---,r JV )= <f>«(h,---, WMn-LJ 

x a(r 2 - 1 2 ) • ■ -a(r N - l N ). (8) 

Here 1 labels the sites of the lattice, and the coeffi- 
cients 0a (li, • ■ • , ljv) are taken to be completely sym- 
metric or antisymmetric. Assuming that we are at the 
minimum of a band with an isotropic effective mass m*, 
we see that the envelope function 0„ (ri, • ■ • , rjv) satis- 
fies a Schrodinger equation with a Hamiltonian given by 
(^) apart from the replacement of m in Hq with m* (the 
mass in Hi is not changed). This is the standard effective 
mass approximation, and it leads to 
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F = y(e z xV)| dzjdl- (V-V>eff( r , r ')L 



where 

Peff(r,r') =iVV/ Q / d 3 r 2 d 3 r 3 • • • (iV 



(9) 



x 4>* a {r,r 2 , - ■ ■ ,r N ) 4> a (r' ,r 2 , - ■ ■ ,r N ) (10) 

is a one-particle density matrix constructed from the en- 
velope functions. 

Unfortunately, the integrand in (|^) is not proportional 
to the actual canonical momentum density j(r), the lat- 
ter having fluctuations on the scale of the lattice con- 
stant b. However, it is possible to prove that when any 
local single-particle operator like the current density is 
averaged over a length scale larger than the localization 
length of the Wannier functions (assumed to be of the 
order of b), but smaller than the characteristic length 
scale over which the envelope functions vary, the expec- 
tation value is correctly given by p c ff(r,r') p3|. In other 



words, the structure on the scale of the lattice constant 
is not described correctly by the envelope functions, but 
coarse-grained quantities are. Now, because the line in- 
tegral of the actual canonical momentum is quantized, it 
is possible to write it as 

Idl j = - / da- [ dzidl-j, (11) 

J a a 2 - ai J ai L JO J a 



where both radii a\ and a 2 are much larger than A and 
their difference is larger than b. The integrals over a and 
z on the right-hand-side of (|ll|) have the effect of averag- 
ing the azimuthal component of j. Therefore, the actual 
circulation is correctly given by the line integral of the 
coarse-grained momentum, which, in turn, is correctly 
given by the envelope functions. Hence, (J9j) leads to the 
transverse force (||), as stated. 

The total transverse force per unit length on the vor- 
tex may be obtained by adding the conventional Lorentz 
force term fL = — - n v c x leading to 



f=pKx(V-v c ) = ^ J d 2 rn(V- v c ) x B. 



Alternatively, we can write this as 



J d 2 r n (v p - v c ) x B + - J d 2 r n (V - v p ) x B, 



(12) 



(13) 



where v p is the velocity of the positively charged substrate, usually taken to be at rest. The first term in (13) is 



the Lorentz force, given by the interaction of the Galilean-invariant current with the magnetic field, while the second 
term is a Magnus force that acts on the substrate. This interpretation is in agreement with the early picture of flux 
line motion in type-II superconductors given by Nozieres and Vincn in the late sixties jl4| : The Magnus force on 
the vortex may be thought as the Kutta-Joukowski hydrodynamic lift force due to the circulation of the electron 
fluid around the vortex. Far from the vortex the circulation is reduced but an increase in the Lorentz force exactly 
compensates this deficiency. The Magnus force reaction is eventually carried away by the positive substrate. 
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APPENDIX A: MEISSNER EFFECT IN THE RPA 
APPROXIMATION 

As mentioned above, it is important to establish that 
the Darwin Hamiltonian provides an accurate model of a 
charged superfluid or superconductor; in particular, that 
it exhibits a Meissner effect with the conventional Lon- 
don penetration depth. 

This can be demonstrated by treating the current- 
current interactions in an RPA-type approximation and 
calculating the current induced by a weak applied vector 
potential, 



^(q)=X u '(q)^ x t(q)- (Al) 

The zero-frequency linear response function x %J is the 
sum of a retarded current-current correlation function 
for a system described by the Hamiltonian Hq + Hi , 
and a diamagnetic term. Our RPA approximation corre- 
sponds to a summation of the diagrams shown in Fig. |l| 
for the imaginary-time current-current correlation func- 
tion n y (iw„, q), which are the most divergent terms as 
q — » 0. In this approximation we find 

X = Xo + XoTx, (A2) 

where Xo ^ s the corresponding response function for the 
system without current-current interactions, as described 
by H , and T^'(q) = 47r(^'|q|- 2 - q l q 3 H^)/c is the 
Fourier transform of T y '(r). 

The response function x y can be used to relate the 
total vector potential A to t = A oxt + A; n( j, the sum of an 
external and induced part, to A cxt itself, 
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^ ot (q) = [1 - T(q) X o(q)]^ cxt (q). (A3) 

Because u> = here, ( |A3[ ) tells us about magnetic screen- 
ing. Now, Xo 1S known to have some general properties 
in the q — > limit, reflecting the presence of off-diagonal 
long-range order. In particular, 

Xo J '(q) = — (^-?V/|q| 2 ) (A4) 
mc 

in this limit [p^|, resulting in 

Atot(q) = (l + Aext(q), (A5) 

where A = (mc 2 /47rn s e 2 )5 is the usual London penetra- 
tion depth. Thus, the current-current interaction term 
leads to a conventional Meissner effect. 
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FIG. 1. RPA approximation for the correlation function 1T-' (iui n , q). Here the solid lines denote the exact Green's functions 
for Ho and the zig-zag lines represent the current-current interaction T lJ (q). 
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